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Abstract 

We discuss applications of the JIMWLK renormalization group equation to multigluon correlations in high energy 
collisions. This includes recent progress in computing the energy dependence of higher point Wilson line correlators 
from the JIMWLK renormalization group equation. We find that the large N c approximation used so far in the 
phenomenological literature is not very accurate. On the other hand a Gaussian finite N c approximation is surprisingly 
close to the full result. We also discuss correlations at large rapidity separations, relevant for the "ridge" correlations 
observed in experiments. 
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1. Introduction 

The physics of high energy hadronic or nuclear collisions is dominated by the gluonic degrees of freedom of the 
colliding particles. These small x gluons form a dense nonlinear system that is, at high enough i/s, best described 
as a classical color field and quantum fluctuations around it. The color glass condensate (CGC, for reviews see [T]) 
is an effective theory developed around this idea. It gives an universal description of the small x degrees of freedom 
that can equally well be applied to small x DIS as to dilute-dense (pA or forward AA) and dense-dense (AA or very 
high energy pp) hadronic collisions. The nonlinear interactions of the small x gluons generate dynamically a new 
transverse momentum scale, the saturation scale Q s , that grows with energy. At high enough energy the color glass 
condensate is thus a one-scale system, characterized by a dominant momentum scale Q s that is hard enough to justify 
a weak coupling calculation. The scale Q s dominates both the gluon spectrum and multigluon correlations. The 
nature of a unique saturation scale as both the typical gluon transverse momentum and as the correlation length l/Q s 
differentiates the CGC qualitatively from the high-x part of the wavefunction. 

The most convenient parametrization of the dominant gauge field is in terms of Wilson lines that describe the 
eikonal propagation of a projectile through it. The Wilson lines are drawn from a probability distribution, whose 
dependence on rapidity is described by the JIMWLK renormalization group equation. It reduces, in a large N c and 
mean field approximation, to the BK [2| equation and further, in the dilute linear regime, to the BFKL one. 



2. Correlations in a dilute-dense collision 

One of the more striking signals of saturation physics at RHIC is seen in the relative azimuthal angle Aip de- 
pendence of the dihadron correlation function, where the Aip as n back-to-back peak is seen to be suppressed in 
dAu-collisions compared to pp collisions at the same kinematics |4|. The CGC description of this correlation starts 
from a large x parton pair propagating eikonally through the target. For a dilute target momentum conservation causes 
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Figure 1 : The JIMWLK result for the quadrupole correlator compared to the Gaussian approximation. Shown are the initial condition (MV model) 
aty = and the result after 5.18 units of evolution in rapidity, for the "line" (left) and "square" (right) coordinate configurations. Figures from 
Ref. (3). 



a peak in the correlation function at A<p w n, whereas for a dense one the intrinsic transverse momentum, of the order 
of 2 S , causes the peak to disappear for semihard momenta. To calculate the matrix element for this process one needs 
target expectation values of products of Wilson line operators, such as the dipole and the quadrupole 

D(x T ,y T ) = ^-Ti-U(x T )UHy T ) Q(x T ,y T ,u T ,v T ) = -^Trf/(x r )f/ t (y r )f/(u r )f/ t (v r ). (1) 

JV C iv c 

For practical phenomenological work it would be extremely convenient to be able to express these higher point cor- 
relators in terms of the dipole, which is straightforward to obtain from the BK equation. In the phenomenological 
literature so far [5| this has been done using a "naive large N c " approximation where the quadrupole is assumed to 
be simply a product of two dipoles. A more elaborate scheme would be a "Gaussian" approximation ("Gaussian 
truncation" in O), where one assumes the relation between the higher point functions and the dipole to be the same 
as in the (Gaussian) MV model. The expectation value of the quadrupole operator in the MV model has been derived 
e.g. in Ref. Q. 

In Ref. [3] the validity of these approximations was studied by comparing them numerically to the solution of the 
JIMWLK equation. As studying the full 8-dimensional phase space for the quadrupole operator would be cumber- 
some, we have concentrated on two special coordinate configurations. The "line" configuration is defined by taking 
Uj = x T and V7- = yj, with r — \x T -y T \ and the "square" by taking x T , y T , Uj-, V7- as the corners of a square with side 
r. 

Our most important results in Ref. [3| for the quandrupole expectation value are shown in Figs. [T] and [2j with a 
comparison of the initial and evolved (for 5.18 units in y) results to the approximations. The MV-model initial con- 
dition v = satisfies the Gaussian approximation by construction. Figure [T] shows that the Gaussian approximation 
is still surprisingly well conserved by the evolution. A possible explanation for this based on the structure of the 
JIMWLK equation has recently been proposed in JHJ . The naive large N c approximation, on the other hand, fails 
already at the initial condition, as shown in Fig. [2] This stresses the importance of the various SU(3) group struc- 
ture constraints violated in this approach. Crucially for the phenomenological consequences, even the characteristic 
length/momentum scale differs by factor ~ 2 from the actual result. 

This result does not yet fully address the effect on the measurable dihadron cross section. For that one must 
convolute the Wilson line operators with the q — > qg splitting wavefunction. This nontrivial numerical task is still 
work in progress, discussed in this conference in Q- 

3. Unequal rapidity correlations and the ridge 

Calculating multigluon correlations in a collision of two dense gluonic systems requires including the JIMWLK 
evolution for both of the colliding projectiles (for a more formal discussion see [ 10 1). Of particular phenomenological 
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Figure 2: The JIMWLK result for the quadrupole correlator compared to the "naive large N c " approximation. Shown are the initial condition (MV 
model) at y = and the result after 5.18 units of evolution in rapidity, for the "line" (left) and "square" (right) coordinate configurations. Figures 
fromRef. (3). 



interest are correlations at large rapidity separations, because they are directly responsible for the "ridge" long range 
rapidity correlations observed in nucleus -nucleus and high multiplicity proton-proton collisions. In spite of this, most 
practical work on the JIMWLK equation has so far concentrated on correlations of Wilson lines at one rapidity. 

The existing applications of the CGC framework to calculations of the ridge IfTTl [121 have used a simplified 
argument based on the MV model, where the change in the effective color charge in one infinitesimal step in rapidity 
is independent of the color charge. This leads to an unequal rapidity correlation function of Wilson line operators 
which is independent of the rapidity separation between the produced particles. Consider for example the following 
correlation function, which is similar to the one appearing in the calculation of the ridge correlation in the fcr-factorized 
approximation IfTTl : 

(D(k r ) v D(k r ) v+Av ) 



(N c 2 - 1) 



1 



(2) 



(£>(k T )y){£>(k T ) y+Ay ) 

where D(kj) is the Fourier-transform of the dipole operator|l}. In the MV model calculation of IfTTl this quantity 
would be one independently of Ay. This naturally leads to a very long range correlation in rapidity between produced 
gluons, as seen in IfTTl fT2l . Genuine JIMWLK evolution should be expected to decorrelate the Wilson line operators, 
at a characteristic rapidity scale Ay ~ l/a s . In Fig.|3]we show preliminary results for the correlation function Q. On 
the left is plotted the correlation function itself at different rapidity separations Ay. There is only a mild dependence 
on the transverse momentum kj, as expected from the MV model calculation. The decorrelation seems, however, to 
be relatively fast. This is further quantified on the right, where we show the result for the decorrelation speed f in an 
exponential decay fit 

/ d 2 k r k 4 r [(D(k r ) v D(k r ) v+Av ) - (£(kr) y ) (£(kr) y+ A y )l 

—. ' '' N - 11 ~ expl-^y). (3) 

/ d 2 k r 4D(k r ) v 

Here the weighting with k 2 is chosen to match the unintegrated gluon distribution ~ k^D(k T ) appearing in a kj- 
factorized calculation of gluon production. In Fig. [3] the decorrelation speed g, explicitly ~ a s , is compared to the 
natural evolution speed A, defined by A = d In Ql/dy, measured in the same JIMWLK simulation. The decorrelation is 
seen to be surprisingly fast. What is more worrying, there seems to be a slight logarithmic dependence on the infrared 
cutoff given by the system size L, which persists over a variety of different initial conditions, lattice sizes, values at 
which the running coupling is frozen in the infrared and other parameters varied in the calculations. At the moment 
there appears to be no natural interpretation of this dependence, leaving the implications for ridge phenomenology 
uncertain. 
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Figure 3: Left: Unequal rapidity correlation J5J as a function of kj / Q s for different rapidity separations. Right: Dependence of the decorrelation 
speed f on Q S L, where L is the linear size of the system. 

4. Conclusion 

We have here argued that multiparticle correlations provide unprecedented experimental insight into the details of 
nonlinear QCD dynamics at small x. In particular, they make it necessary to go beyond the mean-field BK equation 
and use the full JIMWLK equation. Azimuthal angle correlations of particles produced in a collision with a dilute 
probe and dense, saturated target are sensitive to multipoint functions of Wilson lines in the target wavefunction. We 
have shown that a "naive large N" approximation used in the literature, where the quadrupole operator is assumed 
to be a simple product of dipoles, is far from the true finite N Q result. A Gaussian approximation, based on the MV 
model, turns out to be surprisingly close to the result from JIMWLK evolution. We have also discussed preliminary 
results on unequal rapidity correlations which are needed for a proper CGC calculation of the "ridge" correlation in 
high energy collisions. We have shown that at least some observables display a problematic infrared behavior of the 
decorrelation speed in rapidity. 
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